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Lg  Attenuation  and  Source  Studies  using  1982  Miramichi  Data 


by  T.-C.  Shin  and  R,  B.  Herrmann 


Abstract 

Using  data  from  earthquakes  in  the  1982  Miramichi  earthquake  source  zone,  spectral  excita¬ 
tion  and  attenuation  of  the  Lg  phase  is  studied.  With  data  in  the  distance  range  of  135  to  994 
km,  interpretation  is  complicated  by  the  presence  of  high  frequency  Sn  and  Pn  waves  which  inter¬ 
fere  with  the  Lg  wave.  At  the  larger  distances,  the  signal  at  frequencies  above  7  Hz  is  completely 
dominated  by  the  non-Lg  arrivals.  A  frequency  dependent  Lg-Q  is  determined  which  rises  from 
300  at  0.5  Hz  to  about  1400  at  10  Hz  at  which  it  flattens  at  higher  frequency.  The  Sn  coda 
apparent  Q  rises  above  3000  at  frequencies  higher  than  10  Hz.  Seismic  moment  and  corner  fre¬ 
quency  estimates  are  made  using  Lg-Q  corrected  spectra.  The  moment  estimates  compare  well 
with  'hose  obtained  from  long  period  surface  waves  and  short  distance  spectral  estimates.  The  Lg 
corne  •  frequency  estimates  are  substantially  lower  than  the  short  distance  estimates.  This 
discrepancy  is  the  subject  of  discussion,  but  the  Lg  moment-corner  frequency  estimates  do  model 
observed  data  well  when  using  a  Brune  (1970)  source  model  and  the  derived  attenuation  relation. 

INTRODUCTION 

The  1982  sequence  of  earthquakes  near  Miramichi,  New  Brunswick  has  been  the  subject  of 
many  studies  (Choy  et  al.,  1983;  Cranswick  et  al.,  1982;  Mueller  and  Cranswick,  1985;  Saikia  and 
Herrmann,  1985;  Wetmiller  et  al.,  1984)  which  discuss  earthquake  locations,  focal  mechanism  and 
source  dynamics  based  on  teleseismic  and  local  observations.  Regional  data  are  also  available  from 
st.  tions  of  the  Eastern  Canada  Telemetered  Network  (ECTN)  at  distances  of  135  to  994  km  from 
the  focal  region. 

The  ECTN  is  a  digitally  telemetered  network  consisting  of  vertical  component  instruments 
with  a  system  velocity  sensitivity  that  is  essentially  flat  between  corner  frequencies  of  1.0  and 
16.0  Hz  (Hasegawa,  1985).  Data  from  the  ECTN  have  previously  been  used  to  estimate  seismic 
moments  and  corner  frequencies  of  eastern  Canada  earthquakes  (Hasegawa,  1983)  and  to  study 
Lg-wave  attenuation  (Hasegawa,  1985),  Hasegawa  (1985)  used  earthquakes  with  travel  paths  to 
ECTN  stations  that  were  within  or  along  the  edge  of  the  Canadian  shield.  Data  paths  for  sources 
in  the  east  in  the  structurally  complex  Appalachian  province  were  excluded  from  his  study.  The 
purpose  of  this  present  study  is  to  use  ECTN  data  to  examine  the  Lg-wave  attenuation  for  those 
paths  excluded  by  Hasegawa  (1985)  and  then  to  compare  results. 


SPATIAL  ATTENUATION 


Figure  1  is  a  map  of  eastern  Canada  showing  the  location  of  the  Miramichi  source  region 
and  the  sites  of  the  ECTN  stations  used  in  the  study.  The  ECTN  data  are  acquired  at  60  samples 
a  second.  Since  this  study  is  interested  in  the  frequency  dependent  excitation  and  attenuation  of 
the  Lg  signal,  the  observed  signals  are  filtered  through  fifteen  different  Butterworth  bandpass 
filters  whose  high-  and  low-frequency  responses  fall  and  rise,  respectively,  at  24  db/'octave.  The 
filter  center  frequencies  were  in  the  range  from  0.5  to  15,0  Hz  and  the  filter  bandwidth  was  taken 
to  be  0.7  times  the  center  frequency. 


The  Lg  spatial  data  are  fit  by  a  model  of  the  form 

A(r)  =  Ao  r-n  e~v  (D 

where  n  is  a  coefficient  that  is  taken  to  be  n  =  5/6  for  Lg  time  domain  measurements  and 
n  =  1/2  for  frequency  domain  measurements.  The  parameter  q  is  related  to  Q  by 

-  _  ttf  (o\ 

QU 

where  f  is  the  frequency  and  U  is  the  group  velocity  of  the  wave  propagation,  taken  to  be  3  5 
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km  si’c 

To  process  the  data,  different  approaches  were  used  for  time  and  frequency  domain  measure¬ 
ments.  The  time  domain  study  is  performed  by  applying  a  recursive  digital  Butterworth  filter  to 
the  data,  and  peak  amplitudes  are  read  in  the  time  window  corresponding  to  the  Lg  arrival.  To 
make  frequency  domain  measurements,  a  17.07  second  time  window,  corresponding  to  1024  time 
domain  samples,  is  used  about  the  expected  Lg  arrival.  The  signal  is  Fourier  transformed,  the 
instrument  response  removed,  and  the  amplitude  spectra  at  the  desired  center  frequency  is 
obtained  by  averaging  the  amplitude  spectra  between  the  two  corner  frequencies  of  filter. 

To  test  the  methodology  for  interpretation,  we  synthesized  the  vertical  component  Lg-wave 
at  each  seismograph  station  by  adding  up  all  higher  mode  surface  waves  for  the  earth  model  of 
Table  1.  A  source  depth  of  5  km,  a  focal  mechanism  with  60°  dip,  100°  slip  and  200°  strike  were 
assumed.  Synthetic  seismograms  were  generated  at  a  sampling  interval  of  60  samples  per  second. 
Three  sets  of  synthetics  were  generated  corresponding  to  source  corner  frequencies  of  1,  2  and  3 
Hz.  The  signals  included  the  ECTN  instrument  response  in  order  to  make  the  synthetics  compar¬ 
able  to  the  actual  observed  ECTN  data.  The  same  processing  steps  are  then  applied  to  the  data 
sets.  Approximately  700  seconds  of  time  history  were  synthesized  at  each  station. 

The  first  test  was  to  measure  the  peak  Lg  amplitude  as  a  function  of  distance  for  each 
earthquake  using  the  synthetics  generated  for  a  perfectly  elastic  model.  Equation  (1)  is  then  used 
with  y  =  0  to  determine  the  coefficient  n.  This  exercise  led  to  an  n  =  0.834,  which  is  in  agree¬ 
ment  with  previous  assumptions  (Nuttli,  1973)  about  the  spatial  decay  of  the  broad  band  Lg 
wave.  The  next  step  involved  the  use  of  the  bandpass  filtered  signals.  Figure  2  presents  the 
geometrical  spreading  factors  n  obtained  from  the  smoothed  spectral  estimates  and  the  filtered 
time  domain  signals.  The  solid  horizontal  lines  indicate  the  values  of  n  =  5/6  for  the  time 
domain  and  n  =  1/2  for  the  frequency  domain  studies.  The  errors  arise  from  the  fact  that  Lg 
wave  amplitudes  do  not  vary  in  a  uniform  manner  with  distance  and  because  radiation  pattern 
effects  for  the  focal  mechanism  used  can  cause  azimuthal  amplitude  variation  not  accounted  for 
by  the  model.  From  this  exercise,  we  conclude  that  the  bandpass  filters  used  are  not  really  so  nar¬ 
row  that  time  domain  geometrical  spreading  cannot  be  used  for  the  recursively  filtered  data.  This 
is  really  expected  from  surface  wave  theory  because  the  additional  time  domain  decay  is  really 
due  to  the  increasing  duration  of  the  signal  with  distance,  and  this  will  also  be  apparent  in  the 
signals  filtered  as  described  here. 

The  geometrical  spreading  of  the  filtered  time  domain  data  is  not  exactly  the  n  =  5^6 
expected  for  time  domain  data.  However,  considering  the  error  bars,  we  will  use  the  n  =  5/6  for 
the  filtered  time  domain  signal.  If  the  actual  geometrical  spreading  is  somewhat  less,  this  will 
mean  that  the  inferred  7  value  is  somewhat  low  and  the  corresponding  Q  value  somewhat  high. 

Table  2  gives  the  epicentral  coordinates,  identification  code,  and  ECTN  published  magni¬ 
tude  for  each  event  studied.  The  range  of  data  covers  over  three  units  in  short  period  magnitude, 
making  the  data  set  an  excellent  one  for  analysis. 

Prior  to  performing  any  detailed  analysis  it  is  essential  that  the  data  be  examined.  Data  at 
short  distances  were  as  expected,  with  the  Lg  arrival  being  the  largest  signal  in  all  frequency 
bands.  However,  we  observed  an  interesting  phenomena  at  distances  greater  than  500  km.  Figure 
3  illustrates  an  this  feature  of  the  data  for  the  recording  of  event  No.  5  at  the  station  WHO,  an 
epicentral  distance  of  702  km  from  the  source.  The  top  trace  gives  the  broadband  data,  while  the 
other  traces  are  narrow  band  pass  filtered  at  the  indicated  center  frequencies  The  peak  Lg 
arrival  is  indicated  by  the  arrow  at  the  top  of  the  figure.  The  Lg  arrival  is  distinct  at  the  low  fre¬ 
quencies,  but  quickly  disappears  within  the  mantle  Sn  coda,  which  is  itself  riding  on  top  of  the 
mantle  Pn  coda  The  implication  of  this  observation,  is  that  we  can  not  arbitrarily  window  the 
trace  about  the  Lg  arrival  and  expect  the  spectrum  to  tell  us  everything  about  the  Lg  at  all  fre¬ 
quencies  With  this  particular  signal,  we  can  only  use  frequency  up  to  7  0  Hz  if  we  are  able  to 
determine  how  much  the  Lg  rises  above  the  coda  of  the  Sn. 

Equation  (1)  is  easily  linearized  and  adapted  to  multiple  events  The  regression  analysis 
\  i e  1  d s  a  single  \  alue  of  7,  which  is  related  to  Q  by  (2),  and  an  A0  term  for  each  event  proportional 
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to  the  source  excitation  at  each  filter  center  frequency.  The  inferred  Q  values  from  the  time 
domain  and  frequency  domain  measurements  using  all  data  at  all  distances  are  shown  in  Figure  7. 
The  Q  estimates  made  in  each  domain  agree  with  each  other.  The  dependence  of  Lg-Q  with  fre¬ 
quency  can  be  approximately  modeled  by  the  relation  Q(f)  =  (5Q0-550)f°'(6.  An  examination  of 
the  data  indicates  a  distinct  change  in  the  frequency  dependence  of  Q  at  about  7  Hz.  Above  this 
frequency  there  is  a  definite  linear  trend,  while  below  this  frequency  there  is  a  suggestion  that  the 
Lg  Q  is  attempting  to  flatten  at  higher  frequencies. 

Because  of  the  observation  that  the  Lg  arrival  time  window  is  contaminated  by  the  Sn  and 
Pn  coda  at  higher  frequencies  and  large  distances,  the  Q  values  given  above  are  suspect  above  7 
Hz.  To  resolve  this  mixture  of  time  domain  data  types,  two  approaches  were  taken.  The  first  used 
the  data  at  all  distances  for  frequencies  less  than  7  Hz,  and  data  acquired  at  distances  less  than 
600  km  for  the  higher  frequencies.  These  Q  values  are  shown  as  the  squares  for  low  frequencies 
and  as  the  asterisks  at  higher  frequencies  in  Figure  5.  The  circles  in  Figure  5  are  the  time  domain 
results  shown  in  Figure  4a.  The  asterisks  have  large  95%  confidence  limits  because  of  the  reduced 
distance  coverage.  The  asterisks  support  the  flattening  of  Q(f)  at  higher  frequency  suggested  by 
the  low  frequency  data. 

The  second  approach  attempts  to  use  the  Lg  data  at  larger  distances  by  estimating  the  Lg 
amplitude.  To  do  this  we  consider  the  envelope  of  the  Lg  arrival  to  be  sitting  upon  a  predictable 
envelope  of  the  Sn  arrival  and  its  coda.  We  fit  the  Sn  coda  shape  using  an  Aki  and  Chouet  (1975) 
body-wave  scattering  model.  Subtracting  the  estimated  coda  envelope  from  the  observed 
envelope  permits  an  estimate  of  the  peak  Lg  amplitude.  The  modified  data  set  is  used  in  the 
regression  analysis  to  yield  the  Q  estimates  given  by  the  diamonds  in  Figure  5.  This  approach 
yields  Q  values  at  higher  frequencies  intermediate  between  the  other  two  estimates. 

SPECTRAL  ATTENUATION 

Anderson  and  Hough  (1984)  presented  an  interesting  discussion  of  modeling  the  high  fre¬ 
quency  falloff  of  acceleration  spectra  in  which  they  showed  that  a  linear  trend  was  observable  if  a 
semi-logarithmic  plot  of  acceleration  spectra  as  a  function  of  frequency  is  made.  This  linear  trend 
indicates  that  at  high  frequencies  the  acceleration  spectra  can  be  modeled  by  an  exponential 
decay  operator.  They  further  found  that  the  exponential  decay  increased  with  distance  which 
could  be  explained  by  a  distance  independent  site  effect  and  a  distance  dependent  propagation 
effect. 


The  acceleration  spectra  A(f)  can  be  parametrically  described  by  combining  a  Brune  (1970) 
source  spectrum  and  an  attenuation  operator  as 


A(f)  = 


^  c-*ft/Q. 

1  +  (fc/f)2 


(3) 


where  ft  is  the  corner  frequency,  t  is  the  travel  time  of  the  arrival,  and  Q  is  the  apparent  propaga¬ 
tion  Q.  If  the  model  is  correct,  then  a  semi-logarithmic  plot  of  A(f)  vs  f  will  exhibit  a  linear  trend 
at  frequencies  greater  than  the  corner  frequency.  If  the  model  is  incorrect,  the  data  may  exhibit 
an  apparent  linear  trend  if  only  a  narrow  range  of  frequencies  is  considered.  One  way  to  test  the 
appropriateness  of  the  Brune  (1970)  model  is  to  form  the  spectral  ratio  of  two  events  at  the  same 
distance.  The  common  attenuation  effect  divides  out,  and  the  resulting  ratio  will  exhibit  a  f~2 
trend  between  the  comer  frequencies  if  the  event  with  the  higher  corner  frequency  has  a  lower 
displacement  spectral  level  at  high  frequencies,  which  is  usually  the  case.  The  F2  trend  is  a  direct 
consequence  of  the  model  Figure  6  compares  the  displacement  spectral  ratios  of  events  1  to  3  at 
EBN.  135  km,  and  LPQ,  257  km.  from  the  source  zone.  The  solid  line  indicates  an  T*  trend  and 
the  arrows  indicate  the  corner  frequencies  estimated  by  fitting  the  Ao  source  terms  of  the  regres¬ 
sion  analysis  to  a  Brune  (1970)  spectrum.  The  vertical  axis  is  plotted  in  arbitrary  units  and  no 
correction  has  be*m  made  for  noise.  The  noise  level  for  the  smaller  earthquake  is  equivalent  to  the 
signal  level  at  frequencies  less  than  0.5  Hz.  This  exercise  indicates  that  the  omega-squared  spectral 
model  is  reasonable,  although  the  plots  are  not  overwhelmingly  convincing. 

Figure  7  shows  the  results  of  processing  the  acceleration  spectra  of  events  3  and  5  at  the 
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station  GNT.  The  straight  line  indicates  the  inferred  linear  treaa.  Based  on  the  appearance  of  the 
linear  trends,  weights  were  assigned  to  each  slope  determination  as  were  windows  between  which 
the  slopes  were  to  be  determined.  In  this  figure,  data  were  used  in  the  3  -  8  Hz  range,  and  the 
slope  was  determined  to  be  -0.0930  logio  units  per  Hz.  If  we  assume  a  group  velocity  of  3.5  km,  s 
to  the  station,  the  inferred  Q  is  1853  using  data  from  6  events.  Figure  8  presents  the  spectral  Q 
values  determined  from  data  at  each  station.  The  form  of  presentation  points  out  the  spectral  Q 
determined  seems  to  be  very  stable  for  a  wide  range  of  epicentral  distance.  Note  that  the  Q 
values  given  here  are  not  directly  related  to  the  spatial  Q  values  given  in  Figures  4  and  5  since  a 
single  frequency  independent  Q  value  is  used  to  model  the  spectra.  Since  this  Q  is  based  on  the 
slope  which  is  in  turn  dependent  upon  the  high  frequency  spectral  level,  the  spectral  Q  deter¬ 
mined  from  the  slope  is  more  and  indication  of  the  high  frequency  Q  rather  than  a  lower  fre¬ 
quency  Q.  We  note  that  the  spectral  constant  Q  of  about  1900  is  within  the  range  of  possible  Q 
values  near  10  Hz  in  Figure  5. 


CORNER  FREQUENCIES  AND  SEISMIC  MOMENT 

The  ground  vertical  component  Lg-wave  displacement  spectra,  D(f),  determined  by  correct¬ 
ing  the  observed  spectra  by  the  displacement  sensitivity  of  the  instrument,  at  a  distance  r,  is 
modeled  by  a  relation 


D(f)=  M°  M(  In  )i/2e-»ft/Q 
4rrp^3  r0  r 


(5) 


where  M0  is  the  seismic  moment,  p  is  the  density,  3  is  the  crustal  shear  wave  velocity,  r0  is  the 
distance  marking  a  change  in  signal  character  from  a  single  body-wave  arrival  to  surface-wave 
arrivals  beyond  this  point,  and  S(f)  is  a  normalized  Brune  displacement  source  spectrum,  given  by 

s(f)  =  (Mf/g2)-1 

This  relation  was  given  by  Herrmann  and  Kijko  (1983a)  to  predict  the  log-mean  spectral  level. 
We  use  the  Hasegawa  (1983)  values  of  r0  =  100  km,  p  —  2.8  gm/cm3,  and  3  =  3.8  km/sec,  to 
apply  this  relation  to  eastern  Canada.  Since  we  can  take  Q(f)  =  550f°9f>,  approximately,  a  non¬ 
linear  inversion  technique  can  be  used  to  use  (5)  to  determine  M„  and  fc.  Table  3  presents  these 
estimates  of  these  parameters  using  all  available  data.  The  table  entries  give  the  event  number, 
the  catalog  magnitude,  Mm,  the  average  fc  and  its  standard  deviation,  the  log-mean  estimate  of 
the  seismic  moment,  obtained  by  averaging  the  log  moments  and  then  taking  the  antilog,  and  the 
multiplicative  confidence  factor  EM0  to  which  the  moments  are  known.  The  last  column,  N„  gives 
the  number  of  spectra  used  in  the  estimation. 

In  addition  to  listing  the  parameters  fc  and  Mq  which  describe  the  observed  data,  estimates 
of  source  radius,  r0  and  stress  drop,  Act,  are  given  using  the  Brune  (1970)  relations 

_  2.34,3 


and 


Act  — 


7M0 

16r03 


W  e  also  derived  an  ra^  magnitude  by  low  pass  filtering  the  ECTN  trace  to  simulate  a 
seismometer-galvanometer  system  with  a  1.0  Hz  seismometer  and  a  galvanometer  with  0.8  critical 
damping  and  a  1  5  Hz  natural  frequency.  In  order  to  truly  model  a  WWSSN  short  period  system, 
another  lowpass  filter  at  3.0  Hz  would  be  required  to  model  the  effect  of  the  Benioff  seismometer 
inductance,  the  equivalent  response  of  the  simulated  seismograms  is  sufficiently  similar  to  that  of 
the  WWSSN  short  period  instrument  to  be  used  in  determining  the  mLg  magnitude.  We  use  trie 
mlg  definition  of  Herrmann  and  Nuttli  (1982) 

mLg(f)  =  2  94^0.8331ogio(r/ 10)  e  0.4342~/r  —  !og1(>A(r,f).  (6) 


where  r  is  the  epicentral  distance  in  km.  A  is  the  reduced  ground  amplitude  in  microns  at 


I 
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frequency  f,  and  7  is  the  regional  spatial  Lg  attenuation  coefficient  and  here  equals  0.00183  km'1 
at  1.5  Hz.  Herrmann  and  Kijko  (1983b)  showed  that  this  relation  can  be  used  with  instruments 
other  than  the  WNVSSN  short  period  if  used  as  described.  The  values  of  the  computed  mug  are 
also  given  in  Table  3,  and  can  be  compared  to  the  ECTN  catalog  values.  The  computed  mLg 
values  have  standard  deviations  varying  between  0.16  and  0.22  magnitude  units. 

An  independent  check  of  the  derived  sesimic  moments  was  made  by  comparing  the  source 
terms  determined  from  the  regression  of  the  frequency  domain  Lg  spectra  against  distance.  There 
was  good  agreement.  We  also  corrected  each  spectra  for  Q(f),  equalized  for  surface  wave  geometr¬ 
ical  spreading,  and  determined  an  RMS  average  amplitude  spectrum.  These  agreed  very  well  with 
the  spectra  predicted  for  the  source  parameters  of  Table  3. 

Figure  9  compares  the  seismic  moment  -  corner  frequency  values  determined  in  this  study 
with  a  scaling  law  proposed  by  Nuttli  (1983)  and  with  observed  scaling  on  the  Canadian  Shield 
that  Hasegawa  (1983)  obtained  using  ECTN  data.  Figure  10  plots  the  seismic  moment  versus 
magnitude  values.  The  circles  use  the  ECTN  catalog  magnitudes  and  the  asterisks  use  the  mLg 
values  determined  in  this  study.  The  Nuttli  (1983)  proposed  scaling  law  and  the  empirical 
Hasegawa  (1983)  scaling  laws  are  shown  for  reference. 

DISCUSSION 

This  paper  has  presented  a  number  of  source  parameter  and  attenuation  estimates  that 
must  be  evaluated  in  light  of  other  studies.  The  seismic  moment  estimates  are  quite  acceptable. 
Nguyen  (1985)  determined  a  seismic  moment  of  5.9  x  1023  dyne-cm  for  event  1  by  performing  a 
simultaneous  inversion  for  focal  mechanism,  focal  depth  and  seismic  moment  from  long  period 
Love-  and  Rayleigh-wave  data  recorded  in  North  America.  This  agrees  well  with  our  Lg  estimate 
of  5.6  x  1023  dyne-cm.  Saikia  and  Herrmann  (1985)  modeled  three  component  time  history  data 
from  data  within  10  km  of  the  source  of  event  3,  and  estimated  the  seismic  moment  to  be  in  the 
range  of  3.3  -  4.1  x  1020  dyne-cm,  to  which  our  estimate  5.9  x  1020  dyne-cm  agrees  well.  Some  of 
the  events  of  this  study  also  provided  locally  recorded  data  analyzed  by  Mueller  and  Cranswick 
(1985;.  These  were  events  3,  11,  12,  13  14,  15  and  16.  Mueller  and  Cranswick  (1985)  used  spectral 
techniques  to  estimate  the  seismic  moments,  corner  frequencies  and  stress  drops.  Keeping  in  mind 
that  our  seismic  moments  for  events  11  through  16  were  based  on  only  one  vertical  component 
spectra,  the  agreement  of  the  estimates  is  quite  good.  Our  estimates  are  usually  less  than  a  factor 
of  2  greater  than  theirs,  the  only  exception  being  event  15  for  which  our  moment  is  approximately 
2.5  to  3  times  greater  than  theirs.  This  comparison  of  Lg  seismic  moment  estimates  with  other 
independent  estimates  points  out  the  stability  of  Lg  seismic  moment  estimates  when  care  is  taken 
to  account  for  anelastic  attenuation. 

The  value  of  using  Lg-wave  corner  frequency  estimates  to  infer  source  properties  is  currently 
the  subject  of  considerable  discussion  (Mueller  and  Cranswick,  1985).  There  seems  to  be  an 
apparent  contradiction  in  that  the  Lg  inferred  stress  drops  are  usually  lower  than  the  comparable 
estimates  using  short  distance  data.  We  would  like  to  avoid  this  issue  for  the  present  by  stating 
that  the  combination  of  a  Brune  (1970)  source  spectrum  model  and  an  attenuation  operator  suc¬ 
cessfully  fit  ali  observed  spectra  (Appendix  B,  Shin,  1985).  Using  short  distance  data  for  event  3, 
Mueller  and  Cranswick  (1985)  found  S-wave  corner  frequencies  in  the  range  of  6  -  15  Hz,  while 
Saikia  and  Herrmann  (1985)  found  values  in  the  range  of  15  -  20  Hz  for  the  same  data  set.  The 
Lg-comer  frequency  of  Table  3  is  constrained  to  be  near  3.5  Hz.  One  major  difference  in  the  data 
sets  is  that  the  short  distance  spectra  could  not  be  fit  by  a  Brune  (1970)  omega-squared  spectrum, 
but  rather  required  high  frequency  displacement  spectrum  falloff  with  frequency  greater  than  F", 
sometimes  greater  than  F3!  To  further  complicate  the  comparison  of  short  distance  spectral  esti¬ 
mates  to  large  distance  Lg  spectral  estimates,  Cranswick  el  al  (1985)  discussed  the  problems  with 
near  surface  site  resonances,  w  hich  can  either  increase  or  decrease  the  apparent  corner  frequency 
depending  upon  the  resonance  frequency.  One  can  argue  that  the  Lg  may  be  the  more  robust  indi¬ 
cator  of  source  properties  since  it  samples  more  of  the  focal  sphere,  it  has  different  angles  of 
incidence  at  the  receiver  which  may  smooth  out  the  resonance  effects  seen  in  a  single  incident  ray, 
and  by  the  fact  that  its  randomness  of  arrivals  avoids  the  coherence  required  for  resonance  The 
Lg  study  does  support  a  higher  frequency  content  in  the  source  than  does  the  Nuttli  (1983)  scaling 
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model,  as  seen  in  Figure  9. 

The  attenuation  estimates  are  not  without  some  ambiguity  either.  We  have  found  that  it  is 
difficult  to  see  the  Lg  wave  at  high  frequencies  at  large  distances  for  these  travel  paths  crossing 
into  the  Canadian  Shield  because  it  attenuates  faster  than  the  coda  following  the  Sn  and  Pn  (Fig¬ 
ure  3).  This  observation  complicates  the  analysis  of  Lg-Q  but  also  indicates  that  high  frequency 
waves  can  propagate  large  distances  in  the  upper  mantle.  This  observation  is  only  now  being 
made  because  of  seismographs  recording  in  this  frequency  range.  Our  estimate  of  Lg-Q  for  tne 
region  is  given  in  Figure  5.  Until  further  data  are  analyzed  at  high  frequency,  our  preferred 
model  would  be  to  have  the  Lg-Q  rise  sharply  at  low  frequencies  and  then  flatten  to  about  1500- 
2000  at  high  frequencies.  This  would  explain  the  Lg  spectral  Q  of  1800  or  so  observed  at  short  dis¬ 
tances,  and  would  also  admit  the  possibility  of  the  Sn  coda  overwhelming  the  Lg  at  larger  dis¬ 
tances.  The  mantle  Sn  and  its  coda  can  be  presumed  to  be  a  problem  only  beyond  the  distance  at 
which  Sn  first  appears.  Hasegawa  (1985)  modeled  vertical  component  acceleration  data  derived 
from  the  ECTN  recordings  and  concluded  that  the  Lg-Q  ws^  Q(f)  =  900f°2.  This  estimate,  based 
on  a  larger  data  set,  but  a  different  set  of  propagation  paths,  is  within  the  error  bars  of  our  esti¬ 
mates  at  frequencies  greater  than  2.0  Hz.  The  Hasegawa  (1985)  has  better  statistical  properties 
because  of  the  wider  range  of  epicentral  distances  available,  whereas  this  study  used  16  earth¬ 
quakes  recorded  at  16  stations,  with  the  source  receiver  distances  fixed. 

The  seismograms  we  examined  indicated  a  relationship  between  arrivals  expected  from  wave 
theory  and  those  due  to  a  scattering  mechanism.  Following  Dainty  (1981)  the  frequency  depen¬ 
dence  of  the  observed  Lg-Q,  Qob<(f),  can  be  expressed  as 

Qobi(f)‘l  =  Q*b.'1  -r  (Qo..cJ)-’ 

which  states  that  the  observed  Q(f)  is  a  superposition  of  a  frequency  independent  intrinsic, 
absorptive  Q,  and  a  frequency  dependent  scattering  Q.  For  our  data  set  we  would  estimate  the 
intrinsic  Q  to  be  about  1500-2000  and  the  Q0  to  be  about  600. 
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Figure  Captions 

Fig.  1.  Distribution  of  ECTN  stations  and  location  of  the  Miramichi  source  zone  (star). 

Fig.  2.  Comparison  of  Lg  geometrical  spreading  made  by  regressing  peak  amplitudes  of  the  syn¬ 
thetic  data  set  with  distance  using  time  domain,  top,  and  smoothed  spectral  amplitudes  in 
the  frequency  domain,  bottom,  as  a  function  of  the  Butterworth  filter  center  frequency. 

Fig.  3.  Digital  data  from  the  station  WBO  for  event  No.  5.  The  broadband  signal  is  given  in  the 
top  trace,  while  the  filtered  traces  are  shown  below. 

Fig.  4.  Lg-Q  values  determined  from  time  domain  (a)  and  frequency  domain  (b)  analysis.  95  % 
confidence  limits  are  given. 

Fig.  5.  Lg-Q  values  from  time  domain.  Circles  represent  use  of  data  at  all  distances,  asterisks  the 
use  of  data  at  distances  less  than  600  km,  and  the  diamonds  the  use  of  data  corrected  for 
the  Sn  coda. 

Fig.  6.  Displacement  spectral  ratio  of  event  No.  1  to  event  No.  3  recorded  at  stations  EBN  (a)  and 
LPQ  (b).  The  arrows  indicate  the  corner  frequencies  and  the  solid  line  an  P2  trend.  The 
vertical  scale  is  in  arbitrary  units. 

Fig  7.  Example  of  a  semilog  plot  of  acceleration  spectra  at  station  GNT  for  events  No.  3  and  5. 
The  solid  line  indicates  the  mean  slope  corresponding  to  Q=1853  for  this  station.  The 
slope  is  determined  using  data  in  the  3.0  -  8.0  Hz  range. 

Fig.  8.  Map  showing  the  source  region,  star,  and  the  high  frequency  Q  determined  for  each  sta¬ 
tion,  which  is  plotted  adjacent  to  the  station  location. 

Fig.  9.  Plot  of  Lg  comer  frequency,  fc,  versus  seismic  moment,  M,,,  for  the  events  studied.  The 
scaling  laws  proposed  by  Nuttli  (1983)  and  Hasegawa  (1983)  are  shown  for  comparison. 

Fig  10.  Plot  of  seismic  moment,  Mq,  versus  mLf.  The  circles  are  the  ECTN  catalog  values  and  the 
asterisks  are  those  determined  in  this  study.  In  addition  the  Nuttli  (1983)  scaling  law  for 
mid-plate  earthquakes  and  the  Hasegawa  (1983)  relation  for  eastern  Canada  are  given. 
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Mn  :  Lg  magnitude  used  in  ECTN 
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Fig.  2.  Comparison  of  Lg  geometrical  spreading  made  by  regressing  peak  amplitudes 
of  the  synthetic  data  set  with  distance  using  time  domain,  top,  and  smoothed 
spectral  amplitudes  in  the  frequency  domain,  bottom,  as  a  function  of  the 
Butterworth  filter  center  frequency. 
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Fig  7  Example  of  a  semilog  plot  of  acceleration  spectra  at  station  GNT  for  events 
No  3  and  5.  The  solid  line  indicates  the  mean  slope  corresponding  to  Q=I853 
for  this  station  The  slope  is  determined  using  data  in  the  3.0  -  8.0  Hi  range. 
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Fig.  8.  Map  showing  the  source  region,  star,  and  the  high  frequency  Q  determined  for 
each  station,  which  is  plotted  adjacent  to  the  station  location. 


shown  for  comparison, 


Fig.  10.  Plot  of  seismic  moment,  Mq,  versus  m^.  The  circles  are  the  ECTN  catalog 
values  and  the  asterisks  are  those  determined  in  this  study.  In  addition  the 
Nuttli  (1983)  scaling  law  for  mid-plate  earthquakes  and  the  Hasegawa  (1983) 
relation  for  eastern  Canada  are  given. 


